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We study quark confinement physics using lattice QCD. In the maximally abelian 
(MA) gauge, the off-diagonal gluon amplitude is strongly suppressed, and then the 
off-diagonal gluon phase shows strong randomness, which leads to a large effective 
off-diagonal gluon mass, M of j ~ f ,2GeV. Due to the large off-diagonal gluon mass 
in the MA gauge, low-energy QCD is abelianized like nonabelian Higgs theories. 
In the MA gauge, there appears a macroscopic network of the monopole world-line 
covering the whole system. We extract and analyze the dual gluon field £? M from 
the monopole-current system in the MA gauge, and evaluate the dual gluon mass 
as mg = 0.4 ~ 0.5GeV in the infrared region, which is a lattice-QCD evidence 
of the dual Higgs mechanism by monopole condensation. Even without explicit 
use of gauge fixing, we can define the maximal abelian projection by introducing a 
"gluonic Higgs field" <j>(x), whose hedgehog singularities lead to monopoles. From 
infrared abelian dominance and infrared monopole condensation, infrared QCD is 
describable with the dual Ginzburg-Landau theory. In relation to the color-flux- 
tube picture for baryons, we study the three-quark (3Q) ground-state potential 
V3Q in SU(3) C lattice QCD at the quenched level, with the smearing technique for 
enhancement of the ground-state component. With accuracy better than a few %, 
V3Q is well described by a sum of the two-body Coulomb part and the three-body 
linear confinement part a"3Q-L m ini where L m i n denotes the minimal value of the 
total length of the color flux tube linking the three quarks. Comparing with the 
Q-Q potential, we find a universal feature of the string tension as o^q ~ °"qq an d 
the OGE result for the Coulomb coefficient as A$q ~ j^QQ- 

1. Nonabelian Feature in QCD 

Nowadays, quantum chromodynamics (QCD) is established as the funda- 
mental theory of the strong interaction. In spite of the simple form of QCD 
lagrangian 



it is still hard to understand nonperturbative QCD (NP-QCD) phenomena, 
such as color confinement and dynamical chiral-symmetry breaking, in the 
infrared strong-coupling region. In fact, QCD tells us the "rule" of the 
elementary interaction between quarks and gluons, but to solve QCD is 
another difficult problem. 




(1) 
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In this point, QCD may resemble the Rubik Cube, where the rule is 
quite simple, but to solve this puzzle is rather difficult. Like QCD, the 
difficulty of the Rubik cube comes from the non-commutable procedures 
of the "nonabelian" rotational process as shown in Fig.l. Here, a kind of 
"local procedure" on 3 x 3 x 3 subcubes leads to extremely large number of 
configurations, which makes this puzzle difficult and interesting. Then, the 
Rubik cube may be regarded as a miniature of QCD. 1 The main difficulties 



Step 1 : Rotation on 
(vTjj upper side 



Step 1 



Step 2 I 



Step 2: Rotation on 
right side 



\ 



Step 2 



I Step 1 



Figure 1. The difficulty and the interest of the Rubik cube originate from non- 
commutable operations based on the nonabelian feature of the rotational group. The 
configuration after Step 1 and Step 2 depends on the order of these operations. 



of QCD originate from the nonabelian and the strong-coupling features in 
the infrared region below 1 GeV. In the ultraviolet region, the QCD coupling 
is weak and then we can use the perturbative QCD, where the three and 
the four gluon vertices stemming from the nonabelian nature of QCD are 
treated as perturbative interactions. In the infrared region, however, the 
nonabelian and the strong-coupling features are significant. 

In the electro-magnetism, the superposition of solutions is possible, be- 
cause of the linearity of the field equation d^F^ — j v , so that we can 
consider the partial electro-magnetic field formed by each charge and the 
total electro-magnetic field can be obtained by summing up the individ- 
ual field solution. On the other hand, the QCD field equation becomes 
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nonlinear as 

+ig[A^G^\ =f, (2) 

due to the nonabelian feature, and then it is difficult to solve it even at the 
classical level, because the superposition of solutions is impossible. So, we 
cannot divide the color-electromagnetic field into each part formed by each 
quark. Instead, the QCD system is to be analyzed as a whole system. In 
this way, for the analysis of QCD, the nonabelian feature provides one of 
serious difficulties. 

2. Abelianization of QCD in MA Gauge 

The nonabelian nature of QCD in the infrared region is, however, removed 
in the MA gauge. In fact, QCD is reduced into an abelian gauge theory 
with color-magnetic monopoles, keeping essence of infrared nonperturbative 
features. 

2.1. Dual Superconductor Theory and QCD 

In QCD, to understand the confinement mechanism is one of the most 
difficult problems remaining in the particle physics. As is indicated by 
hadron Regge trajectories and lattice QCD calculations, the confinement 
phenomenon is characterized by one- dimensional squeezing of the color- 
electric flux and the string tension a ~ lGeV/fm, which is the key quantity 
of confinement. 

On the confinement mechanism, Nambu first proposed the dual super- 
conductor theory for quark confinement, 2 based on the electro- magnetic 
duality in 1974. In this theory, there occurs the one-dimensional squeezing 
of the color-electric flux between quark and anti-quark by the dual Meiss- 
ner effect due to condensation of bosonic color-magnetic monopoles. How- 
ever, there are two large gaps between QCD and the dual superconductor 
theory. 3-5 

(1) The dual superconductor theory is based on the abelian gauge the- 
ory subject to the Maxwell-type equations, where electro-magnetic 
duality is manifest, while QCD is a nonabelian gauge theory. 

(2) The dual superconductor theory requires condensation of color- 
magnetic monopoles as the key concept, while QCD does not have 
color-magnetic monopoles as the elementary degrees of freedom. 

These gaps may be simultaneously fulfilled by taking MA gauge fixing, 
which reduces QCD to an abelian gauge theory including color-magnetic 
monopoles. 
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2.2. MA Gauge Fixing and Relevant Mode for Confinement 

In Euclidean QCD, the maximally abelian (MA) gauge is denned so as to 
minimize the total amount of the off-diagonal gluons, 3 " 5 

R oS [A»(-)]= [ d 4 xtr{[D„,H][D„,H?} = e ^ f d 4 x^K(x)\ 2 , (3) 

by the SU(iV c ) gauge transformation. Here, we have used the Cartan de- 
composition, A^x) = A^(x) ■ H + J2 a A ^( x ) Ea - Since thc SU(iV c ) co- 
variant derivative operator = + ieA^ obeys the adjoint gauge trans- 
formation, the local form of the MA gauge condition is easily derived as 
[H, tf]]]=0. 3 - 5 

In the MA gauge, the gauge symmetry G = SU(A r c ) loca i is reduced into 
H = U(l)^°~ 1 1 x Weylfj° bal , where the global Weyl symmetry is a subgroup 
of SU(iV c ) relating the permutation of N c bases in the fundamental repre- 
sentation. In thc MA gauge, off-diagonal gluons behave as charged matter 
fields like Wjf- in the Standard Model, and provide the color-electric current 
in terms of the residual abelian gauge symmetry. In addition, according to 
the reduction of the gauge symmetry, color-magnetic monopoles appear as 
topological objects reflecting the nontrivial homotopy group 3-8 

n 2 (su(7v c )/u(i) Ar - 1 ) = iMua)^- 1 ) = z£-\ (4) 

in a similar manner to similarly in the GUT monopole. Here, the global 
Weyl symmetry and color-magnetic monopoles are relics of nonabelian na- 
ture of QCD. 

Thus, in the MA gauge, QCD is reduced into an abelian gauge theory 
including color-magnetic monopoles, which is expected to provide a theo- 
retical basis of the dual superconductor theory for quark confinement. Fur- 
thermore, recent lattice QCD studies show remarkable features of abelian 
dominance and monopole dominance for NP-QCD in thc MA gauge. 

(1) Without gauge fixing, all the gluon components equally contribute 
to NP-QCD, and it is difficult to extract relevant degrees of freedom 
for NP-QCD. 

(2) In the MA gauge, QCD is reduced into an abelian gauge theory in- 
cluding the electric current and the magnetic current which 
forms the global network of the monopole world-line covering the 
whole system. In the MA gauge, lattice QCD shows abelian domi- 
nance for NP-QCD (confinement, chiral symmetry breaking 9 , gluon 
propagators 10 ): only the diagonal gluon, which remains an abelian 
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gauge field, is relevant for NP-QCD, while off-diagonal gluons do 
not contribute to NP-QCD. 
(3) By the Hodge decomposition, the diagonal gluon is decomposed 
into the "photon part" and the "monopole part" , corresponding 
to the separation of j M and k^. In the MA gauge, lattice QCD 
shows monopole dominance 911 ^ 12 for NP-QCD: the monopole part 
(kfj, ^ 0, jfj, = 0) leads to NP-QCD, while the photon part 
(iju 7^ 0, /c M = 0) seems trivial like QED and does not contribute 
to NP-QCD. For example, on the Q-Q potential, the purely lin- 
ear confinement potential appears in the monopole part, while the 
Coulomb potential appears in the photon part like QED. 13 

In fact, by taking the MA gauge, the relevant collective mode for NP-QCD 
can be extracted as the color-magnetic monopole. 4 ' 5 

3. Strong Random Phase of Off-diagonal Gluon in MAQCD 

To find out essence of the MA gauge, we study the feature of the off-diagonal 
gluon field A± = -±= (A* ± iAfJ in the MA gauge in SU(2) lattice QCD. 

In SU(2) lattice QCD, the SU(2) link variable is factorized as U^(s) = 
M^(s)u )1 (s), according to the Cartan decomposition SU(2)/U(1) 3 x U(l) 3 . 
Here, u^(s) = exp{ir 3 # 3 (s)} 6 U(l) 3 denotes the abelian link variable, and 
the abelian projection is defined by the replacement as U^(s) — > u M (s). The 
off-diagonal matrix M M (s) S SU(2)/U(1) 3 is parameterized as 

M ^ s >~ e ~ [ie^hme^s) cos^(s) ) ' (b) 

In the continuum limit, Xn( s ) coincides with the off-diagonal gluon phase 
as A^(x) = e ±%Xf> ( x )\Ap(x)\, and ^(s) is proportional to the off-diagonal 
gluon amplitude as O^s) = -^a\eA^{s)\ with the lattice spacing a. In the 
MA gauge, the diagonal element cos^(s) in M^(s) is maximized by the 
SU(2) gauge transformation, e.g. (cos ^(s)) ma — 0.93 at [3 = 2.4. Ac- 
cordingly, the off-diagonal gluon A^(x) = e ±lx ^^\A^(x)\ has two relevant 
features in the MA gauge. 3-5 

(1) The off-diagonal gluon amplitude \A^(x)\ (or | sin# p (s)| on lattices) 
is strongly suppressed by SU(iV c ) gauge transformation in the MA 
gauge. 

(2) The off-diagonal gluon phase Xn( x ) tends to be random, because 
Xfj,( x ) is n °t constrained by MA gauge fixing at all, and only the 
constraint from the QCD action is weak due to a small accompany- 
ing factor \A^\ or | sinfyj. 
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Figure 2. (a) The distribution P(Ax) of the difference Ax = |Xm( s ) — X^( s + ^)l (mod7r) 
in the MA gauge plus U(l)3 Landau gauge at /3=0 (a = oo, thin line), (3=1.0 (a ~ 0.57fm, 
dotted curve), 0=2.4 (a ~ 0.127fm, solid curve), ,8=3.0 (a ~ 0.04fm, dashed curve), (b) 
The off-diagonal gluon contribution to the Wilson loop Wq 1 v.s. the perimeter length 
L = 2(R + T) in the MA gauge on 16 4 lattice with /3=2.4. The thick line denotes the 
theoretical estimation in Eq.(10) with the microscopic input (ln{cos 6> M (s)})ma — —0.082 
at = 2.4. 

Now, we consider the behavior of Ax = |Xm( s ) ~ Xp( s + ^')l( nl0 d7r) in the 
MA gauge. If the off-diagonal gluon phase Xn( x ) 1S a continuum variable, as 
lattice spacing a goes to 0, Ax — a\d u x^\ must go to zero, and hence P(A%) 
approaches to the <5-function with the peak at A^ — 0. However, as shown 
in Fig. 2(a), P(A%) is almost a-independent and almost flat. These features 
indicate the strong randomness of the off- diagonal gluon phase Xni x ) m the 
MA gauge. Then, Xp.( x ) is approximately a random angle variable in the 
MA gauge. a 

3.1. Randomness of Off-diagonal Gluon Phase and Abelian 
Dominance 

Within the random-variable approximation 3 " 5 for the off-diagonal gluon 
phase X/i( s ) m the MA gauge, we analytically prove abelian dominance of 
the string tension. Here, we use 

(e* (s) )MA * ±- r dxv(s) exp{i X/t (*)} = 0. (6) 
2tt Jo 

In the Wilson loop (W C [U}) = (trLTcU^s)) = (trn c {M M (s)u M (s)}}, the 
off-diagonal matrix M At (s) is reduced to a c- number factor, 

Af M (a) -»cosfl„(«) 1, (7) 



a However, as shown in Sect. 6, near the monopole, there remains a large amplitude 
of |sin# M (s)| even in the MA gauge, and Xm( s ) * s strongly constrained so as to make 
the QCD action small. Hence, Xm( s ) cannot be regarded as a random variable near 
monopoles. 
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and then the SU(2) link variable U^(s) is reduced to a diagonal matrix as 

U„(s) = M„(s)«„(s) -► cos^(s)u M (a), (8) 

after the integration over Xm( s )- For ^ ne R x T rectangular C, the Wilson 
loop Wc[f ] in the MA gauge is approximated as 

(W C [U}) = (trnf =1 {M w ( Sj )u w ( Si )}) * (n?=i cose w («i) •trn j i =1 w, 1 .( Sj )}MA 
~ (exp{Ef =1 ln(cos^ i (s i ))})MA (W c [«])ma (9) 
with perimeter L = 2(i? + T) and the abelian Wilson loop 
W c [u] = trll^^^^Sj). Replacing J2i=i^ n { cos (°^( s i))} b y its average 
L(ln(cos^ A1 ))MA in a statistical sense, we derive the perimeter law of the 
off- diagonal gluon contribution to the Wilson loop as 

Wf = (W c [U])/(Wc[u})ma * exp{i(ln(cos^)) M A}, (10) 

which is confirmed in lattice QCD, as shown in Fig. 2(b). Near the con- 
tinuum limit, we find also the relation between the macroscopic quan- 
tity Wq S and the microscopic quantity of the off-diagonal gluon amplitude 

(|eA±| 2 ) M A as 

Wf ee (Wc[U])/(Wc[u}) MA ~ exp{-i^(| e A±| 2 } M A}. (11) 

In this way, perfect abelian dominance for the string tension, C7gu(2) = OAbei, 
is analytically derived within the random-variable approximation for the 
off-diagonal gluon phase in the MA gauge. 

3.2. Strongly Random Phase and Large Mass of 
Off-diagonal Gluons 

As another remarkable fact, strong randomness of off-diagonal gluon phases 
leads to rapid reduction of off- diagonal gluon correlations. In fact, if X/j-( x ) 1S 
a complete random phase, Euclidean off-diagonal gluon propagators exhibit 
^-functional reduction, 

(A+(x)A-(y)) MA = {\A+(x)\\A- u (y)W^-^) M K 

= (\A±(x)\ 2 ) MA 6 flv 6 4 (x-y), (12) 

which means the infinitely large mass of off-diagonal gluons. Of course, the 
real off-diagonal gluon phases are not complete but approximate random 
phases. Then, off-diagonal gluon mass would be large but finite. In this 
way, strong randomness of off-diagonal gluon phases is expected to provide 
a large effective mass of off-diagonal gluons. 
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4. Large Mass Generation of Off-diagonal Gluons in MA 
Gauge : Essence of Infrared Abelian Dominance 

Using SU(2) lattice QCD, we actually investigate the Euclidean gluon prop- 
agator Gf v {x-y) = (A a ^(x)Al(y)) (a, 6=1,2, 3) and the off-diagonal gluon 
mass M Q ff in the MA gauge. As for the residual U(l)3 gauge symmetry, we 
take U(l)3 Landau gauge, to extract most continuous gluon configuration 
under the MA gauge constraint and to compare with the continuum the- 
ory. The continuum gluon field A®(x) is derived from the link variable as 
U„(s) = exp{iaeA%(s)^}. 

We show in Fig. 3(a) the scalar-type gluon propagators GjL(r) and 
G+"(r) = (A+(x)A-(y)) = ^{G^r) + G%„ (r)}, which depend only on 
the four-dimensional Euclidean distance r = J (x M — y^) 2 , in SU(2) lattice 
QCD with 2.2 < (3 < 2.4 with various sizes (12 3 x 24, 16 4 , 20 4 ). We find 
infrared abelian dominance for the gluon propagator in the MA gauge: only 
the abelian gluon A 3 (x) propagates over the long distance and can influence 
the long-distance physics. 4 ' 5,10 

Since the four-dimensional Euclidean propagator of the massive vector 
boson with the mass M takes a Yukawa-type asymptotic form as 

3 M 1 3M 1 / 2 P ~ Mr 

<Wr) = ^~^Mr) + W2 S\x - y) * (13) 

we investigate the effective mass M D s of off-diagonal gluons A^ (x) from the 
slope of the logarithmic plot of r 3 / 2 G+~(r) ~ exp{— M gr} in Fig. 3(b). 

(1) The off-diagonal gluon A^ (x) behaves as a massive field with a large 
mass about 1 GeV for r > 0.2fm in the MA gauge. 

(2) From the fitting analysis of the lattice QCD data with r > 0.2fm, 
the off-diagonal gluon mass is evaluated as M Q s ~ 1.2 GeV in the 
MA gauge. 

We perform also the mass measurement of off-diagonal gluons from the 
temporal correlation of the zero- momentum projected operator 0^(t), 

r+"(r) = (O+(t)O-(0)), 0±(r) = J dx A±(x,r), (14) 

in lattice QCD in the MA gauge plus U(l) 3 Landau gauge. We find the off- 
diagonal gluon mass M Q g ~ 1.2GeV again from the slope of the logarithmic 
plot of r+-(r) in SU(2) lattice QCD with 2.3 < (3 < 2.35 with 16 3 x 32 
and 12 3 x 24. 
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Thus, off-diagonal gluons acquire a large effective mass M s — 
1.2GeV in the MA gauge, which is essence of infrared abelian 
dominance. 4,5 ' 10 In the MA gauge, due to the large effective mass M Q g ~ 
I.2GeV, off-diagonal gluons can propagate only within a short range 
as r ^ M~ s — 0.2fm, and becomes infrared inactive like weak bosons in 
the Standard Model. Then, in the MA gauge, off-diagonal gluons A^ can- 
not contribute to the infrared NP-QCD, which leads to infrared abelian 
dominance. 4 ' 5 ' 7 ' 10,14 



[fm' 



r M GJf) l ^" 2 ) 12 J X24,16 J 20 4 r * W /r -<0) 



' OIT-diuykniLil 'i\m.i\\ 



I Ji.il'ihi.iI gluon 



" ° 8 '° [fl 




(a) 



(b) 



(c) 



Figure 3. (a) The scalar-type gluon propagator G^ M (r) as the function of the four- 
dimensional distance r in the MA gauge in SU(2) lattice QCD with 2.2 < f3 < 2.4 with 
various sizes (12 3 X 24, 16 4 , 20 4 ). (b) The logarithmic plot of r 3 / 2 G^ M (r) v.s. r. The off- 
diagonal gluon propagator behaves as the Yukawa- type function, Gp^i ~ cxp ^ 3 ^" ff r ^ ■ 

(c) The logarithmic plot of the temporal correlation Via (t) = {O^ (t)O^T (0)) as the 
function of the temporal distance r in SU(2) lattice QCD with 2.3 < /3 < 2.35 with 
16 3 X 32 and 12 3 X 24. From the slope of the dotted lines in (b) and (c), the effective 
mass of the off-diagonal gluon A J is estimated as M g ~ 1.2GeV. 



5. Lattice-QCD Evidence of Monopole Condensation 

In the MA gauge, there appears the global network of the monopole world- 
line covering the whole system as shown in Fig. 4(a), and this monopole- 
current system (the monopole part) holds essence of NP-QCD. Using SU(2) 
lattice QCD, we examine the dual Higgs mechanism by monopole conden- 
sation in the NP-QCD vacuum in the MA gauge. 4 ' 5 

Since QCD is described by the "electric variable" as quarks and gluons, 
the "electric sector" of QCD has been well studied with the Wilson loop 
or the inter-quark potential, however, the "magnetic sector" of QCD is 
hidden and still unclear. To investigate the magnetic sector directly, it is 
useful to introduce the "dual (magnetic) variable" as the dual gluon field 
B^, which is the dual partner of the diagonal gluon and directly couples 
with the magnetic current 
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Owing to the absence of the electric current in the monopole part, 
the dual gluon can be introduced as the regular field satisfying (d A 
B)nv = *Fp U and the dual Bianchi identity, A B)^ v = j v = 0. By 

taking the dual Landau gauge d^B^ = 0, the field equation is simplified 
as d 2 B fl — d a *F a ^ — fc M , and therefore we obtain the dual gluon field B^ 
from the monopole current fc M as 

= (d-%)( X ) = ~J ^(^2 ■ (15) 

Since the dual gluon B^ is massive under monopole condensation, we inves- 
tigate the dual gluon mass ms as the evidence of the dual Higgs mechanism. 

First, we put test magnetic charges in the monopole-current system 
in the MA gauge in SU(2) lattice QCD, and measure the inter-monopole 
potential Vjf(r) to get information about monopole condensation. Since the 
dual Higgs mechanism provides the screening effect on the magnetic flux, 
Vm (r) is expected to be short-range Yukawa- type, if monopole condensation 
occurs. The potential between the monopole and the anti-monopole can be 
derived as 

V M (R) = - lim hn(W D (R,T)), (16) 

T — *oo 1 




Inter-Monopole Distance r [frn] r E [im] 



(a) (b) (c) 

Figure 4. The SU(2) lattice-QCD results in the MA gauge, (a) The monopole world- 
line projected into R 3 on the 16 3 X 4 lattice with (5 = 2.2 (the confinement phase). There 
appears a global network of monopole currents covering the whole system, (b) The intcr- 
monopolc potential Vm(i") v.s. the 3-dimensional distance r in the monopole-current sys- 
tem on the 20 4 lattice. The solid curve denotes the Yukawa potential with tub = 0.5GeV. 
The dotted curve denotes the Yukawa-type potential including the monopolc-sizc effect, 
(c) The scalar-type dual-gluon correlation ln(r|/ 2 (i3 (J (a;)B fl (j/)}MA) as the function of 
the 4-dimensional Euclidean distance te on the 24 4 lattice. The slope corresponds to 
the dual gluon mass nig. 
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using the dual Wilson loop Wo as the loop-integral of the dual gluon, 4,5 

W D {C) = exp{i| £ dx^} = cxp{z^ J j do-^F^}, (17) 

which is the dual version of the abelian Wilson loop 

W Ahcl (C) = exp{i| £ dx^} = exp{i| J J da^F^} (18) 

and the test monopole charge is set to be e/2. 

In Fig. 4(b), we show Vm{t) in the monopole part in the MA gauge. 4 ' 5 
Except for the short distance, the inter-monopole potential is well fitted by 
the Yukawa potential 

V M(r) = J^—, (19) 
Air r 

and thus the magnetic screening is observed. In the MA gauge, the dual 
gluon mass is estimated as ttib — 0.5GeV from the infrared behavior of 
V M (r). 

Second, we investigate also the Euclidean scalar-type dual gluon prop- 
agator (B t ^(x)B fl (y)) ma as shown in Fig. 4(c), and estimate the dual gluon 
mass as tub = 0.4 ~ 0.5 GeV from its long-distance behavior. 5 

From these two tests, the dual gluon mass is evaluated as raj = 0.4 ~ 
0.5 GeV, and this can be regarded as the lattice-QCD evidence for the dual 
Higgs mechanism by monopole condensation at the infrared scale. 4 ' 5 

To conclude, lattice QCD in the MA gauge exhibits infrared abelian 
dominance and infrared monopole condensation, and therefore the dual 
Ginzburg-Landau (DGL) theory 8 ' 17-23 can be constructed as the infrared 
effective theory directly based on QCD in the MA gauge. (See Fig. 5). 



6. Monopole Structure in terms of the Off-diagonal Gluon 

Let us compare the QCD-monopole with the point-like Dirac monopole. 
There is no point-like monopole in QED, because the QED action diverges 
around the point-like monopole. The QCD-monopole also accompanies 
a large abelian action density, however, owing to cancellation with the off- 
diagonal gluon contribution, the total QCD action is kept finite even around 
the QCD-monopole. 3 ^ 5 

To see this, we examine the QCD-monopole structure in the MA gauge 
in terms of the action density using SU(2) lattice QCD. 3 " 5 From the SU(2) 
plaquette P^ {2) {s) and the abelian plaquette P^ e \s) E U(l) 3 C SU(2), 
we define the "SU(2) action density" S^ (2 \s) = 1 - |trP^ (2) (s), the 
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"abelian action density" S^ el (s) = 1 — itrP^, bel (s) and the "off-diagonal 
gluon contribution" S°^,(s) = S^^ 2 \s) — Sff el (s). In the lattice formalism, 
the monopole current k l _ l (s) is defined on the dual link, and there are 6 
plaquettes around the monopole. Then, we consider the local average over 



From Lattice QCD to Dual Ginzburg-Landau Theory 

— Infrared Effective Theory directly based on QCD — 



QCD : SU(3) C Nonabelian Gauge Theory 



Maximally Abelian (MA) Gauge Fixing 

( partial gauge fixing ) [ G.'t Hooft . NPB 190('8 1 )455 ] 



U(l>3 x U(l)g Abelian Gauge Theory+ QCD-monopole 



Lattice QCD studies 

• Abelian Dominance 

Only diagonal gluon is relevant 
forNP-QCD 

• Monopole Condensation 



7I 2 [SU(3)/{U(1) 3 XU(1) 8 }]=Z2 

hedgehog color-magnetic monopole 
configuration 



[cf . GUT - monopole] 



- Dual Ginzburg-Landau Theory 



£dgl= - jtr(^B v - d v B^) 2 +tr[?<? M + X? [&„ + gB p X\ 
-AtrCrtT-* 2 ) 2 

Bp. = Bftl& B*r&, : dual gluon field 

X = XaEa''- QCD-monopole field 
4 jr 

£> = — : dual gauge coupling constant 

e 

X : coupling of monopole self-interaction 

v : imaginary mass of monopole — > monopole condensate 



Dual Gauge Symmetry is spontaneously broken instead of Gauge Symmetry 

Figure 5. Construction of the dual Ginzburg-Landau (DGL) theory from lattice QCD 
in the maximally abelian gauge. 
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the 6 plaquettes around the dual link, 

1 1 

S(s,fi) = -^Yl \ £ ^(3~f\S a p{s + mj). (20) 

a/37 ™=0 

We show in Fig. 6(b) the probability distribution of the action densities 
*Ssu(2) i ^Abei and S g around the QCD-monopolc in the MA gauge. 



Abelian magnetic monopole 




S(iu) ' si.) 

(a) (b) (c) 

Figure 6. (a) The total probability distribution P(S) on the whole lattice and (b) the 
probability distribution Pk(S) around the monopole for SU(2) action density Ssu(2) 
(dashed curve), abelian action density SAbel (solid curve) and off-diagonal gluon con- 
tribution S ff (dotted curve) in the MA gauge at /3 = 2.4 on 16 4 lattice. Around the 
QCD-monopole, large cancellation between SAbcl a,w ^ ^ofi keeps the total QCD-action 
small, (c) The schematic figure for the QCD-monopole structure in the MA gauge. The 
QCD-monopole includes a large amount of off-diagonal gluons around its center as well 
as the diagonal gluon. 

We summarize the results on the QCD-monopole structure as follows. 

(1) Around the QCD-monopole, both the abelian action density SAbel 
and the off-diagonal gluon contribution S s are largely fluctuated, 
and their cancellation keeps the total QCD-action density 5 l gu(2) 
small. 

(2) The QCD-monopole has an intrinsic structure relating to a large 
amount of off-diagonal gluons around its center like the 't Hooft- 
Polyakov monopole. At a large scale, off-diagonal gluons inside 
the QCD-monopole become invisible, and QCD-monopoles can be 
regarded as point-like Dirac monopoles. 

(3) From the concentration of off-diagonal gluons around QCD- 
monopoles in the MA gauge, we can naturally understand the local 
correlation between monopoles and instantons. 3 ~ 5,11,12,15,le In fact, 
instantons tend to appear around the monopole world-line in the 
MA gauge, because instantons need full SU(2) gluon components 
for existence. 
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7. Gluonic Higgs and Gauge Invariant Description of MA 
Projection 

In the MA gauge, the gauge group G = SU(A C ) is partially fixed into its 
subgroup H = U^ 1 x Wcyff}° bal , and then the gauge invariance becomes 
unclear. b 

In this section, we propose a gauge invariant description of the MA 
projection in QCD. Even without explicit use of gauge fixing, we can natu- 
rally define the MA projection by introducing a "gluonic Higgs scalar field" 
4>(x). For a given gluon field configuration {A^(x)}, we define a gluonic 
Higgs scalar 4>(x) — Q(x)HW (x) with Cl(x) G SU(iV c ) so as to minimize 

R[$(-)}= J d^tv^Jix)]^,^}. (21) 

We summarize the features of this description as follows. 3 

(1) The gluonic Higgs scalar <f)(x) does not have amplitude degrees of 
freedom but has only color-direction degrees of freedom, and 4>(x) 
corresponds to a "color-direction" of the nonabelian gauge connec- 
tion Dp averaged over /j, at each space-time x. 

(2) Through the projection along (j>(x), we can extract the abelian 
U(l) Jv c-i su b-gauge-manifold which is most close to the origi- 
nal SU(7V C ) gauge manifold. This projection is manifestly gauge 
invariant, and is mathematically equivalent to the ordinary MA 
projection. 3 

(3) Similar to D^, the gluonic Higgs scalar <fi(x) obeys the adjoint 
gauge transformation, and is diagonalized in the MA gauge. Then, 
monopoles appear at the hedgehog singularities of <j>(x) as shown in 
Fig.7. 3 ' 8 

(4) In this description, infrared abelian dominance is interpreted as 
infrared relevance of the gluon mode along the color-direction 4>{x), 
and QCD seems similar to a nonabelian Higgs theory. 

8. The 3Q Ground-State Potential in SU(3) Lattice QCD 

In relation to the color-flux-tube picture for baryons, 24 ' 25 we study the 
three-quark (3Q) ground-state potential in SU(3) C lattice QCD at the 



b In Refs.[3,12], we show a useful gauge-invariance criterion on the operator Oma defined 
in the MA gauge: //Oma defined in the MA gauge is H -invariant, Oma also invariant 
under the whole gauge transformation of G. 
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Figure 7. The correlation between the gluonic Higgs scalar field <j>(x) = <t> a (x)^- and 
monopoles denoted by dots in SU(2) lattice QCD with /3 = 2.4 and 16 4 . The arrow 
denotes the color direction of ((ft 1 (x) , (j> 2 (x) , <fi 3 (x)) in the SU(2) internal space. The 
monopole in the MA gauge appears at the hedgehog singularity of the gluonic Higgs 
scalar o(x) in the Landau gauge. 

quenched level. 26-28 In contrast with a number of studies on the Q-Q po- 
tential using lattice QCD, there were only a few preliminary lattice-QCD 
works for the 3Q potential. 29-31 

To begin with, let us consider the potential form in the Q-Q and 3Q sys- 
tems with respect to QCD. In the short-distance limit, perturbative QCD 
is applicable and the Coulomb-type potential appears as the one-gluon- 
exchange (OGE) result. In the long-distance limit at the quenched level, 
the flux-tube picture 24 ' 25 would be applicable from the argument of the 
strong-coupling limit of QCD 32 , and hence a linear confinement potential 
proportional to the total flux-tube length is expected to appear. Indeed, 
lattice QCD results for the Q-Q ground-state potential are well described 
by 

^QQM—^+'W + CqQ (22) 

at the quenched level. In fact, Vqq is described by a sum of the short- 
distance OGE result and the long-distance flux-tube result. 

Also for the 3Q ground-state potential V3Q , we try to apply this simple 
picture of the short-distance OGE result plus the long-distance flux-tube 
result. Then, the 3Q potential V3Q is expected to take a form of 

V3Q = — ^3Q ^ 1 — - — I- + <J 3 QL min + C 3 q, (23) 

2<J ' J 1 

with L m i n the minimal value of total length of flux tubes linking the three 
quarks. 

Similar to the derivation of the Q-Q potential from the Wilson loop, 
the 3Q static potential V3Q can be derived from the 3Q Wilson loop W3Q 
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V 3Q = - km i ln(W 3Q ), W 3Q = \,e ahc e a , hld ufufu^ (24) 

1 — >oo 1 o! 

with U k = Pexp{i 5 f T dx^A^x)} (k = 1,2,3) (see Fig.8). Here, the 3Q 
Wilson loop W3Q is defined in a gauge-invariant manner. The 3Q gauge- 
invariant state is generated at t = and is annihilated at t = T. For 
< t < T, the three quarks are spatially fixed in R 3 . 




Figure 8. The 3Q Wilson loop W3Q. The 3Q state is generated at t = and is annihi- 
lated at t = T. The three quarks are spatially fixed in R 3 for < t < T. 

Physically, the true ground state of the 3Q system would be expressed 
by the flux tubes instead of the strings, and the 3Q state which is expressed 
by the three strings generally includes many excited-state components such 
as flux-tube vibrational modes. Since the practical measurement of (W3Q) 
is quite severe for large T in lattice QCD calculations, the smearing tech- 
nique for ground-state enhancement 27,28 is practically indispensable for the 
accurate measurement of the 3Q ground-state potential V3Q . 

The standard smearing for link-variables 33 is expressed as the iteration 
of the replacement of the spatial link-variable Ui(s) (i = 1,2,3) by the 
obscured link-variable Ui(s) G SU(3) C which maximizes 

Re tr(u}(s){aUi(s) + ]T U± j (s)U i (s±j)Ul j (s + i)}'), (25) 

with a real smearing parameter a and J7_^(s) = UUs — /t). The n-th 

smeared link-variables U^ n \s) (n = 1, 2, .., iV smcar ) are iteratively defined 
starting from uj?\s) = U^(s) as 

U^(s) = U^ n - 1 \s) (i= 1,2,3), U { 4 n \ S )^U 4 ( S ). (26) 

This smearing procedure keeps the gauge covariance of the "fat" link- 
variable ul?\s) properly. In fact, the gauge invariance of F(ujT\a)) 
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is ensured if F(C/ M (s)) is a gauge-invariant function. Since the fat link- 
variable Jjjr\s) includes a spatial extension, the "spatial line" expressed 
with £/^(s) physically corresponds to a "flux tube" with a spatial exten- 
sion. Therefore, if a suitable smearing is done, the "spatial line" of the 
fat link-variable is expected to be close to the ground-state flux tube. We 
search reasonable values of the smearing parameter a and the iteration 
number N smeal in the lattice QCD calculation. 



Table 1. The best fit coefficients in Eq.(23) for the 3Q potential and those in Eq.(22) 
for the Q-Q potential in the lattice unit. 







A 


C 


3Q 


0.1524(28) 


0.1331(66) 


0.9182(213) 


QQ 


0.1629(47) 


0.2793(116) 


0.6203(161) 



Table 2. Lattice QCD results for the 3Q potential V^q for 16 patterns of the 3Q system, 
where the three quarks are put on (i, 0, 0), (0, j, 0) and (0, 0, k) in R 3 in the lattice unit. For 
each 3Q configuration, Vjq is measured from the single-exponential fit (W3Q) = Ce~ V 3Q T 
in the range of T listed at the fourth column. The statistical errors listed are estimated 
with the jackknife method, and x 2 /-^DF is listed at the fifth column. The fitting function 
V 3 *q in Eq.(23) with the best fitting parameters in Table 1 is also added. 





V3Q 


C 


fit range 


xV^df 


v 3 % 


" v& 


(0,1,1) 


0.8457( 38) 


0.9338(173) 


5 


-10 


0.062 


0.8524 


-0.0067 


(0, 1,2) 


1.0973( 43) 


0.9295(161) 


4 


- 8 


0.163 


1.1025 


-0.0052 


(0,1,3) 


1.2929( 41) 


0.8987(110) 


3 


- 7 


0.255 


1.2929 


0.0000 


(0,2,2) 


1.3158( 44) 


0.9151(120) 


3 


- 6 


0.053 


1.3270 


-0.0112 


(0,2,3) 


1.5040( 63) 


0.9041(170) 


3 


- 6 


0.123 


1.5076 


-0.0036 


(0,3,3) 


1.6756( 43) 


0.8718( 73) 


2 


- 5 


0.572 


1.6815 


-0.0059 


(1.1,1) 


1.0238( 40) 


0.9345(149) 


4 


- 8 


0.369 


1.0092 


0.0146 


(1,1,2) 


1.2185( 62) 


0.9067(228) 


4 


- 8 


0.352 


1.2151 


0.0034 


(1,1,3) 


1.4161( 49) 


0.9297(135) 


3 


- 7 


0.842 


1.3964 


0.0197 


(1,2,2) 


1.3866( 48) 


0.9012(127) 


3 


- 7 


0.215 


1.3895 


-0.0029 


(1,2,3) 


1.5594( 63) 


0.8880(165) 


3 


- 6 


0.068 


1.5588 


0.0006 


(1,3,3) 


1.7145( 43) 


0.8553( 76) 


2 


- 6 


0.412 


1.7202 


-0.0057 


(2,2,2) 


1.5234( 37) 


0.8925( 65) 


2 


- 5 


0.689 


1.5238 


-0.0004 


(2,2,3) 


1.6750(118) 


0.8627(298) 


3 


- 6 


0.115 


1.6763 


-0.0013 


(2,3,3) 


1.8239( 56) 


0.8443( 90) 


2 


- 5 


0.132 


1.8175 


0.0064 


(3,3,3) 


1.9607( 93) 


0.8197(154) 


2 


- 5 


0.000 


1.9442 


0.0165 



We show lattice QCD results for the 3Q ground-state potential. We 
generate 210 gauge configurations using SU(3) C lattice QCD Monte-Carlo 
simulation with the standard action with (3 — 5.7 and 12 3 x 24 at the 
quenched level. The lattice spacing a ~ 0.19 fm is determined so as to 
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reproduce the string tension as tr=0.89 GeV/fm in the Q-Q potential Vqq. 
Here, the pseudo-heat-bath algorithm is used for updating, and the gauge 
configurations are taken every 500 sweeps after a thermalization of 5000 
sweeps. 

We measure the 3Q ground-state potential V3Q using the smearing tech- 
nique, and compare the lattice data with the theoretical form of Eq.(23). 
Owing to the smearing with a=2.3 and A^ smear =12, the ground-state com- 
ponent is largely enhanced, and therefore the 3Q Wilson loop (W^q) com- 
posed with the smeared link-variable exhibits a single-exponential behavior 
as (Wsq) ~ e~ V3QT even for a small value of T. For each 3Q configuration, 
we extract Vg'g from the least squares fit with the single-exponential form 

(W 3 Q> - Ce- V «* T (27) 

in the range of T m ; n < T < T max listed in Table 2. Here, the fit range of T 
is chosen such that the stability of the "effective mass" 

V(T) = H(Wsq(T))/(W 3Q (T + 1))} (28) 

is observed in the range of T m j n < T < T max — 1. 

For each 3Q configuration, we summarize the lattice QCD data VJq as 
well as the prefactor C in Eq.(27), fit range of T and xV^df in Table 2. 
The statistical error of V^q is estimated with the jackknife method. We 
find a large ground-state overlap as C > 0.8 for all 3Q configurations. 

Now, we consider the potential form of V3Q. We show in Table 1 the 
best fit parameters in Eq.(23) for V3Q. We compare in Table 2 the lattice 
QCD data V$ with the fitting function in Eq.(23) with the best fit 
parameters. The three-quark ground-state potential Vsq is well described 
by Eq.(23) with accuracy better than a few %, although x 2 /A^df = 3.76 
seems relatively large, which may reflect a systematic error on the finite 
lattice spacing. (The fitting with A-type flux-tube ansatz suggested in 
Rcfs. [29,31,34] is rather worse and shows unacceptably large xV-^df = 10.1 
even for the best fit.) 

By comparing the coefficients (<T3q,A3q) with (cqqj^qq) in Table 1, 
we find a universal feature of the string tension, CT3Q ~ 0qq, and the OGE 
result for the Coulomb coefficient, ^4 3 q ~ \ Aqq. 
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